Construction des objets géographiques
Création de la couche géographique des stations
Il s’agit de transformer le tableau de données au format “long” en objet géographique “points” au moyen de ses colonnes de coordonnées (CoordXSiteHydro et CoordYSiteHydro). Le système de projection est Lambert 93.
Vérification
Chargement de la couche shp des SAGEs
Si l’on dispose d’une couche shapefile des SAGEs, il peut être intéressant de les croiser avec les données d’étiage, par exemple pour produire des tableaux de bord par SAGE.
chemin <- "//svfcvin2/DFS/COMMUNS/REGIONS/BRE/DR/dr35_IG_metier/CATALOGUE/EAU/GESTION_TERRITOIRE/SAGE/sage_2016_DIR2_shp.shp"
sages <- lire_couche_sages(fichier_shp = couche_sage, scr = 2154)#> Reading layer `sage_2016_DIR2_shp' from data source `D:\Pascal\boulot\onde\donnees_geographiques_reference\sage_2016_DIR2_shp.shp' using driver `ESRI Shapefile'
#> Simple feature collection with 43 features and 12 fields
#> geometry type: MULTIPOLYGON
#> dimension: XY
#> bbox: xmin: 123905 ymin: 6559835 xmax: 609715.2 ymax: 6892897
#> epsg (SRID): NA
#> proj4string: +proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +units=m +no_defsVérification
Jointure stations - SAGES
Pour associer à chaque station les données sur son SAGE, jointure spatiale :
Comme en Bretagne - Pays de Loire la couche des SAGE est approximative, il y a des pbs à l’issue de la jointure \(\Rightarrow\) un nettoyage de la couche est nécessaire.
On recherche d’éventuels doublons là où les SAGEs se chevauchent. Pour ce faire, on visualise en rouge les points qui sont des doublons (= situés sur + d’un SAGE).
doublons <- stations_onde_geo$CdSiteHydro %>%
table %>%
as.data.frame() %>%
filter(Freq > 1) %>%
.[,1] %>%
as.character()
doublons_geo <- stations_onde_geo %>%
filter(CdSiteHydro %in% doublons)
mapview(stations_onde_geo, viewer.suppress = FALSE, legend = FALSE) +
mapview(sages, zcol = "NOM", alpha = 0.1, legend = FALSE) +
mapview(doublons_geo, color = "red", legend = FALSE)Il apparaît que le seul problème est le Frout (J3104013), mal affecté à cause du chevauchement entre deux périmètres de SAGE. On corrige manuellement. On garde aussi pour visualisation des stations ONDE hors SAGE mais dans les régions de l’étude.
